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Abstract 

We study the stochastic homogenization of the system 



-diva'' = p 

Ph 

where 4>n is a sequence of convex stationary random fields, with p-growth. 

We prove that sequences of solutions (a 11 , vl 1 ) converge to the solutions of a 

deterministic system having the same subdifferential structure. The proof 

relies on Birkhoff 's ergodic theorem, on the maximal monotonicity of the 

subdifferential of a convex function, and on a new idea of scale integration, 

recently introduced by A. Visintin. 
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^ H 1 Introduction and main result 

In this paper we study the behaviour, as 77 — > 0, of the solutions u v '(• , u) , a v '(• , ui) 
/\ of the problem 



— diver'' (x,ui) = P{x) in Q x Vt, 

a v (x,u) e d(f) ri (Vu n (x,uj),x,oj) in Q x Q, (1.1" 

u v (x,u) =0 on dQ x Q, 
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where Q is a bounded domain in IR m and (Q, J 7 , /i) is a probability space. The map 
£ ~~ ► 0(^i w ) is convex for /z-a.e. u£fi and satisfies classical p-growth and coercivity 
conditions for 1 < p < +00. The oscillating function <p v : W 71 x Q x Q — y R is 
defined as 

<j> v (£,x,u) := <j)(£,T x / v u), 

for a dynamical system T^ : Q — >■ $1 on R m . 

A simple example. We refer to Section 2.1 for the general definitions and results 
regarding probability, but in order to have a clearer picture of the setting we 
anticipate a simple one-dimensional example. 

Let p — 2, m — 1, Oi,a2,?i,?2 £ R, Q = (<7i?<?2), and consider a random 
function onMxQ which in every point x & Q has probability P G (0, 1) to have 
the form 0i(£,x) = ai£ 2 and probability 1 — P to have the form (p2{C,,x) = a2^ 2 . 
The idea of homogenization is to approximate this random function by partitioning 
Q into intervals of length 77 > 0, defining the random function 0^ independently 
on each interval, and letting i] tend to zero. Let 

M := {u : K. — > {a±, 02} : 00 is constant on every interval (n, n+1), n G N} , 

then the dynamical system T x can be chosen as the shift operator 



T T : L°°(M) ->■ L 



oc 



- .T 



W(-) i->- w(-+Xj, 



and the set fi is given by the set of all functions obtained from functions in M by 
a shift x G R, i.e. Q = M x T, where T = R/Z is the 1-dimensional torus. The 
probability measure on Q is then the product measure of the measure induced by 
P, 1 — P on M, times the Lebesgue measure on T (for further examples and details 
see, e.g., [14, Section 7.3, "Random Structure of Chess-board Type"] or [7, pages 
349-350]). Finally, for all u which are continuous in 0, define 

<l>T,(t,x,w) :=uj(x/ri)\£\ 2 . 

Rescaling the parameter x of the dynamical system by I/77 corresponds, for 
vanishing 77, to a finer and finer mixing of the realizations of <fi. Under hypothesis 
of ergodicity for T with respect to Q, it is then natural to expect the limit system 
to be deterministic. We prove that for almost every u G Q, sequences (a 71 ,^) of 
solutions of (1.1) converge, as rj — > 0, to the solutions of a subdifferential inclusion 
governed by a deterministic convex function 0o : ^ m — > R- 



Regarding the example above, in the case P = 1/2, it is well known (see, e.g., 
[7, page 350]) that the limit function is given by 

MO = c\0 2 , 

where c = 20102/(01 + 0-2) is the harmonic mean of a\ and a 2 - 
We turn to the main result of the paper. 

Theorem 1.1. Let (Q, J 7 , //) be a probability space with an ergodic dynamical sys- 
tem T x : Q — > Q. Let Q C R m be a bounded domain, letpE (1, +00) , p' = p/{p—l), 
(P)v>o C W^'P'iQ), with p -»■ / strongly in W~ l *\Q). Let c Q ,c l) c 2 > and 
(f) : R m x Q -> R be such that 

u; i->- 0(£, u;) is 7 -measurable V £ 6f ra , (1.2a) 

£ 1— t>0(£,w) is convex for fi-a.e. u £ Q, (l-2b) 

-co + ci|^| p < 0(^,o;) < c + c 2 |C| p /or/i-fl.e.uGO,V(GM m (1.2c) 

Define 

(j) v (£,x,uj) := (f>(Z,T x/rl u), 

then 

(i) for every i] > i/iere exists a couple (u v , a v ) satisfying 

u"(-,w) G ^(Q), ^"(-,w) G L P '(Q) for/2-a.e. u G ft, 

-diva'' = f m W , - 1 *'(Q), /or //-a.e. w G ft, 

o~ v (x, ui) G d(f) v ('Vu 11 (x, w),x, uj) a.e. in Q x ft, (1.3) 

■u ,? (x, w) =0 a.e. on dQ x ft. 

(ii) For fi-a.e. u> G ft the family of solutions (u v (-,u),cr v (-,u)) v> o is weakly 
compact in W^ P {Q) x LP'(Q). 

(Hi) For fi- a.e. u G ft, for every subsequence (^fc)fceN such that 

a vk _! a * y;ea£% in L P '(Q), 
u Vk ^u* weakly inW lj? (Q), 

the limit point (a*,u*) is a solution of 

—diva* = / in Q, 

a* G <90 o (Vw*) in Q, (1.4) 

u* = on dQ. 



The homogenized function 0o : R m — >■ R is given by 

<j>o(Z):=m£\ [<!>(£ + u(u),u) dp: ueL p pot (n), fudfi = o). (1.5) 

(See Section 2.3 for the definition of Lp Qt (fl).) 

1.1 Literature 

We remark that the results of the present paper were already found by Messaoudi 
and Michaille in [17, 18], using the method of epiconvergence of integral function- 
als, in the case of nonconvex integrands with growth p > 1. This result was then 
generalized by Abddaimi, Michaille, and Licht [1], to the case of quasiconvex inte- 
grands with linear growth. As we detail in the next subsection, the main interest 
of the present paper is then to recover the known results with a new, simpler tech- 
nique, which does not involve the convergence properties of integral functionals. In 
this section we review the main literature concerning stochastic homogenization. 

The subject of homogenization has been widely studied since the early works 
of Babuska [2], De Giorgi and Spagnolo [10], Tartar [31, 32], Bensoussan, Lions, 
and Papanicolau [3], Sanchez-Palencia [26]. Restricting ourselves to the case of 
stochastic homogenization (or homogenization in stationary ergodic media), we re- 
call that Kozlov [15] and Papanicolau and Varadhan [22, 23] where the first to 
study second order linear elliptic PDEs. T-convergence techniques for the study 
of random integral functionals were employed by Dal Maso and Modica in [8, 7]. 
In the monograph [14] Jikov, Kozlov, and Oleinik collect results regarding (among 
the other topics) homogenization of linear elliptic problems, elliptic problems in 
perforated random domains and homogenization of random lagrangians with non- 
standard growth conditions via T-convergence, with applications to elasticity. 

More recently, the mathematical analysis community showed a renewed interest 
in stochastic homogenization of different nonlinear problems. As examples we cite 
the papers on parabolic operators [29, 30], the works on Hamilton- Jacobi equations 
[16], on fully nonlinear elliptic equations [6, 5], one-dimensional plasticity [27], and 
on variance estimates for the homogenized coefficients of discrete elliptic equations 
[13]. We also mention a first extension of Nguetseng's two-scale convergence [21] 
to the stochastic setting [4]. 

We note that problem (1.3) was studied in [9] in the setting of periodic homog- 
enization by means of Mosco-convergence and periodic unfolding. In the proof we 
employ a different method of homogenization but we adopt the same variational 
setting and convex analysis tools. 



1.2 Discussion 

As we noted above, the form of the homogenized function <j) is already well-known 
in the literature. The real interest of this paper lies rather in a new proof which 
is direct, in the sense that it does not make use of T-convergence, G-convergence, 
epiconvergence, or two-scale convergence, showing instead simultaneous conver- 
gence of the two components of the solution (o~ v , u v ). The proof is also completely 
self-contained, taken into account that some probability tools like ergodicity (some- 
times traded for independence at large distances, as in [7]) are intrinsecally part 
of the problem. Finally, we apply to the stochastic setting a new idea of homoge- 
nization introduced by Visintin in [34] for the periodic setting (under more general 
hypothesis). This method, named "scale integration," deals with the integration 
with respect to y of the relation 

w(x,y) E dtp(y(x,y),x,y) onQxY, 

where d<p(-, x, y) is a cyclically monotone mapping on a separable real Banach space 
X, Q C M m , and Y is the m- dimensional unit torus. Relying on the splitting of 
the space-regularity that is at the basis of compensated compactness, and on the 
properties of Fenchel-Legendre transform, the method yields 



w(x,y)dy E dip [ v(x,y)dy,xj on Q, 

where <po is the solution of a minimization problem in a suitable subspace V C X 

<Po(Z,x) := inf / ip(£ + u(y), x,y) dy. 



uev j Y 

1.2.1 Description of the proof 

The main ideas underlying the proof of Theorem 1.1 can be summarized as follows 
• We study the auxiliary problem 

— divs = 0, s e <90(Vu>), 

set in the probability space only. Weyl's decomposition of L 2 (Q) spaces into 
potential and solenoidal field provides the correct functional setting for the 
variational formulation and allows for a rigorous definition of the derivatives 
of s and w. Convexity of (f> ensures the existence of a solution. 



By Birkhoff's ergodic Theorem 2.2, for a.e. u G Q the realizations 

s v (x) = s(T x/r) u), Vw v (x) = Vw(T x/v uj) 
converge to their expected values 



E(s) = [ sdfi, E(Vw) = [ V 
Jn Jn 



w dfi. 



• Exploiting the scale integration method, we can relate the expected values 
by the homogenized function <po defined in (1.5) 

E(s) G d<j) (E(Vw)). 

• By compensated compactness (see Lemma 2.2) we can pass to the limit in 
the monotonicity inequality 

(a v - s v , VvP - Vw v ) > 

and recover (1.4) by maximal monotonicity of <90o- 

1.2.2 Further questions 

We remark that in some of the steps above we use tools that allow for more 
general hypothesis than the ones assumed in this work. It would be interesting to 
understand if the scale integration method could be applied to recover the same 
general results which can be obtained via T- or epiconvergence, and if it could yield 
new results. We present a few examples of directions of extension. 

Growth p — 1. All the properties which rely only on convexity and monotonicity 
are still valid, but since L 1 is not reflexive, problems will arise, e.g., from the lack 
of compactness of W 1,1 , which leads to a setting in spaces of bounded variation. 
This problem was solved via epiconvergence, see [1]. 

Nonconvex 0. The scale integration method does not require the integrand <p 
to be convex, so in principle it could be applied also in this case. A number of 
questions arise nonetheless, starting from the lack of lower-semicontinuity for the 
related integral functional. See the results in [17, 18]. 



Time-dependent problems. One of the main physical applications of equations 
(1.3) is to problems of nonlinear viscosity like 

-diva = /, 
d t (V s u) Ed<j)(a). 

Since scale integration applies also to time-dependent problems, we hope to be able 
to use the present result as a starting point for the study of this kind of problems 
in the future. A similar use of an auxiliary problem and oscillating test function 
method was used, in the case of periodic homogenization, in [28]. 

1.2.3 Plan of the paper 

The paper is organized as follows. Section 2 contains all the notation, definitions, 
and the main results regarding ergodicity, Weyl's decomposition, compensated 
compactness, and convex analysis. The proof of Theorem 1.1 is given in Section 3. 

2 Preliminaries 

Notation: for all open sets E C M. m , let T>(E) be the set of all infinitely differen- 
tiable functions with compact support in E, and let T>'(E) be its dual space, that 
is, the distributions on E. 

2.1 Ergodicity 

For this part of classical theory we follow the exposition in [14, Section 7], to which 
we refer for further details and examples. 

Let (Q, J 7 , /i) be a probability space, with a-algebra J 7 and probability measure 
fx. An m- dimensional dynamical system is a family of mappings 

T x : n ->■ il VxGR™, 

satisfying 

1. Group property. Tq — I (where / is the identity), 

T (x+y) =T x T y \fx,yeR m . 

2. Mass invariance. The mappings T x : Q — > Q preserve the measure /jond, 
that is, for every x G M. m and every measurable set E G J 7 we have 

T X E G F, fi{T x E) = fi{E). 



3. Measur ability. For any measurable function / : Q — > W 1 , the function / : 
]R m x Q — > IR n defined as f(x,ou) := f(T x u) is also measurable. 

A random field on W 11 is a mapping £ : IR m x fi — >• IR n which associates to every 
x G lR m a random variable £(x, ■) with values in R n . A random field is said to be 
stationary if for any finite set xi, . . . ,Xk G M. m and any /i G 1R" 1 the distribution of 
the random vector 

t(xi + h),...,£(x k + h) 

does not depend on h. In particular, £ is stationary if it can be represented in the 
form 

£(x,u>) = f(T x u>), (2.6) 

where / is a fixed random variable on Q, and T is a dynamical system. The 
converse statement is analyzed, e.g., in [11]. 

Let a dynamical system T be given, a measurable function / defined on Q is 
said to be invariant if 

f(T x u) = f(u) V(x, W )el m , 

almost surely in Q. A dynamical system T is said to be ergodic if every invariant 
function is constant almost surely in Q. 

Let u I—)- f(oo) be a measurable function on Q. For a fixed u G Q the function 
£ i— > f(T x u) is said to be a realization of f. We often denote the realization of / 
by /(•)■ 

Let p > 1, then we denote by L p (f2) the usual space formed by the equivalence 
classes of J 7 - measurable functions / : Q — > lR m , such that 

\f(uj)\vdn(uj) <+oo, 
<n 

and by L°°(Q) the space of measurable essentially bounded functions. 
Lemma 2.1. The following properties hold: 

• if f G £ p (f2) ; p G [l,+oo[ ; £/ien /i-almost all realizations of f belong to 

• if fn —* f strongly in L P (Q), then there exists a subsequence (/ nfc )fcgN such 
that for p,- a. e. u G Q we have f nk — > f strongly in Lf oc (lR m ) m ; where f nk ( x ) — 
f nk (T x u), f(x) = f(T x u). 



Proof. Using the measure-preserving property of T, for all x G H 

ll/n - /Ilip (n) = J 1/nH " /MI'd/iM 

|/„(T x o;)-/(7;w)|P^(w). 



Let If CC M m , then by the measurability property of T and Fubini's theorem 

HA " /IIW) = nri f f \M T <*>) - f(T x ^dn(uj)dx 

v ' \ K \ JKJn 

= ih I I \fn(T x u)-f(T x u)\ p dxd^u). 
\ K \ JvlJk 

This implies that /i-almost all realizations of / belong to Lj' 0C (R m ) m , and since 
ll/n — f\\ P LP(n) ~ * 0> we can fi nc ^ a subsequence n^ such that 



^n fc H := / \fn k (T x uj) - f(T x u)\dx 
Jk 

converges to zero /i-almost everywhere. □ 

Let / G Ll 0C {R m ) m . A number M{f} is called the mean value of / if 



lim — / f[-)dx = M{f} 

for any measurable bounded set K C M m . Denoting i^ := {x G 1R" 1 : r/x G X}, 
the mean value of / can also be written, after rescaling, as 

lim nT7?T / f(x)dx = M{f}. 

r?->o \K\ J K 

I I «/ x\^ 

The expected value of a random variable / on $1 is defined as 

E(f):= f f{uj)dix. 



JCl 

Theorem 2.2. (Birkhoff's Ergodic Theorem, [14, Lemma 15.1]) Let f G L 1 (Q). 
Then, for a.e. u G $1, the mean value of the realization f(T x u) exists and it satisfies 



M{f(T x u)}dfi = E(f). 
'n 

If the system T is ergodic, then the mean value does not depend on ui almost surely 
and it satisfies 

M{f(T x uu)} = E(f) forfi-a.e.ueQ. 



Remark 1. Birkhoff ' s theorem tells us that almost every realization f v (x) = f(T x / v u) 
satisfies 

lim ,f fr,(x)dx = E(f), 

and it is important to notice that this holds for every measurable bounded set 
K C lR m . This entails in particular that if / G L P (Q), then 



fv - 1 E (f) weakly in Lf c 



ocK w»r. 



2.2 Div-curl Lemma 

Compensated compactness is the fundamental idea for passing to the limit as rj — > 
in a product of weakly converging functions o"' 7 ,^. 

If (j, £ G L 2 (Q) m , and {o - ' 7 }, {V} are two sequences of vector fields in L 2 (Q) m 
such that 

a" -^ a, e" -» e weakly in L 2 {Q) m , 

it is in general false that 

a" ■ e v -» a • e, 

in any sense. Consider, for example, a^x) = ^(x) = sin(x/r)) in L 2 (0, 1). The 
theory of compensated compactness, introduced by L. Tartar and F. Murat in 
[32, 19, 33], gives general conditions on the derivatives of <j v ,e v , and on a function 
F in order to conclude that 

iV,e")-F(<r,e) 

in the sense of distributions. The following statement is one of the consequences 
of compensated compactness theory. 

Lemma 2.2 (Div-Curl lemma, [20]). Let p G (l,+oo) ; p' := p/(p — 1). Let 
ai,ae L p '(Q) m and e r ',ee L p (Q) m such that 

a v ^a weakly in V '{Q) m , e r ' -» e weakly in L p \Q) m . 

In addition, let f G W~ 1,p (Q) and assume that 

curie'' = 0, 

diva* 7 — > / strongly in W~ 1,P '(Q). 

Then 

a v ■ e" A a ■ e inV(Q). 

10 



2.3 Weyl's decomposition 

Weyl's decomposition theorem, or more correctly "Peter-Weyl's" theorem, was 
originally stated for compact topological groups [24]. In the context of stochastic 
homogenization, it is used to provide an orthogonal decomposition of L 2 (Q) into 
functions, the realizations of which are vortex-free and divergence-free, in the sense 
of distributions (see, e.g., [14, Lemma 7.3]). For p G [1, +oo[, Weyl's theorem can be 
generalized to a relation of orthogonality between subspaces of L P (Q) and LP (Q). 
Let v = (v 1 , . . . ,v m ), v l G L^iW 71 ), i — 1, . . . ,m, be a vector field. We say 
that v is potential (or vortex-free) in W 71 if 

v i {x)^^-v j {x)^^dx = \/i,j = l,...,m, V^eV(R m ). 

(2.7) 
Recall that every vortex free vector field in a simply connected domain is a gradient 
of some function u. In IR m it means that there exists u G W l g P {W n ) such that 
v = Vw. We say that a vector field v is solenoidal (or divergence-free) in W 71 if 



Yl v i (x)^r^-dx = V^GP(l m ). (2.8) 

t-^ '™~ ax 4 



i=l 



Now let us consider vector fields on the probability space (Q, J 7 , p). A vector 
field / G L P (Q) is called potential if p- almost all its realizations x \- > f(T x u) 
are potentials, in the sense of (2.7). The space of potential vector fields on Q is 
denoted by lZ ot (Q). A vector field / G L P (Q) is called solenoidal if ^-almost all 
its realizations x i— ¥ f(T x u) are solenoidal, in the sense of (2.8). The space of 
solenoidal vector fields on fi is denoted by L P ol (Q). 

Lemma 2.3. The spaces ZZ ot (fi),l£ ol (fi) satisfy 

Lp 0t (Q) and L p s6[ (Vt) are weakly closed subspaces of L P (Q). (2-9) 

E(a ■ e) = E(a) ■ E(e) for all a G L p ' ol (Q), e G L p pot (Q). (2.10) 

Proof. Let us prove (2.9) for L P ol (Q), the proof for L p t (D,) is identical. Let (u n ) C 
L p ol (Q), and assume that u n — > u, strongly in L P (Q). Then, by definition of L p ol (Q), 
/i-almost every realization u n satisfies 

J2 [ u n {x)^^dx = V V GP(l m ). (2.11) 

. . .IJUrn UXi 



1=1 
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By Lemma 2.1 there exists a subsequence (rik) such that u nk — y u, strongly in 
L p loc {VL). Passing to the limit as k — y oo in (2.11) yields u G Lg 0l (fi), and therefore 
L p ol (Tl) is strongly closed in L P (Q). Being L P ol (Q) a linear subspace of L P (Q), 
strongly closed is equivalent to weakly closed. □ 



Set 



v; = {/a; t (n) : E(f) = o}, 

H P o, = {/eLL(fi) : E(f) = 0}. 



It holds 



L p pot m = V p pot © R m , L P JQ) = V s p ol © R m . 

Lemma 2.4. [14, lemma 15.1] Let p G [1, +oo) ; p' = p/(p — 1). The relations 

(v p p ot ) ± = \& © R m , o&) ± = v;' ot © R m 

hold in the sense of duality between the spaces L P (Q) and L p (Q). 

Remark 2. If v G lZ ot (Q), as noted above, for /z-a.e. u G Q we can find a function 
w such that V x u(x,u) = v(T x u), but this does not imply that u = u(T x ou), i.e., 
that u is a stationary random field. In particular, we are not allowed to apply 
Birkhoff's theorem to the sequence u(x/r),oo). 

2.4 A review of convex analysis 

We review some basic facts of convex analysis. We recall the statements in the 
case when (X, | • \x) is a reflexive Banach space, with dual space (X', | • \x r ) and 
duality pairing "(•, •)", having in mind the applications X = R m , X = L P (Q). 
Let ip : X — y R U {+oo}, the domain of tp is 

dom((p) := {e G X : (p(e) < +00} . 

If p is proper (i.e., dom(<p) 7^ 0), the Legendre-Fenchel transform (or Young- 
Fenchel transform, or conjugate function) </?* is given by 

V?* :X^MU{+oo}, a^sup{(a,e) - (p(e)}. (2.12) 

sex 

The growth of the conjugate function tp* is related to the growth of tp by the 
following result. 
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Lemma 2.5. Let cq,ci, Ci > 0, p G [1, +oo], p' — p/{p — 1). If ip satisfies 

-c + c 1 \e\ p x <ip(e)<c + c 2 \e\ p x Ve G X, 

then there exist ci, C2 > snc/i £/ia£ </?* satisfies 

-c + c^alx, < V * (a) < c + c 2 \cr\ x , V a e X' . (2.13) 

The subdifferential of </?, computed in the point e G dom(ip), is the set 

fy>(e) = { a eX' such that </?(£) > <p(e) + (a,£-e) V(Gl}. 

Some useful properties of convex functions are summarized in the following lemma, 
for a proof we refer to [12]. 

Lemma 2.6. Let ^:X-yRU +oo be proper, then, V e G X , V a G X' ; 

(%) ip* is convex, lower-semicontinuous , and dom(<p*) ^ 0. 

(%%) ip(e) + ip*(a)>(a,e). 

(Hi) a G d<f(s) 4 eG dip* (a). 

(iv) [<p(e) = <p**(e), e G dip* (a)] ^ a e d<p(e). 

(v) e G dip(a) <£> (p(e) + ip*(a) = (a,e). 

The equality in (v) is also known as Fenchel's equality, while (ii) is referred to as 
(Young-)Fenchel's inequality Note that for a proper, convex, lower-semicontinuous 
function <p, implication (iii) becomes an equivalence. 

Define the orthogonal complement V L of a closed subspace l^clas 

V 1 - := {a G X' : Me G V, (a, e) = 0} . 

Theorem 2.3 (Fenchel-Rockafellar). Let V C X be a closed subspace, let ip be 
proper, convex, lower-semicontinuous . Assume that there is a point of V where ip 
is continuous, then 

inf ip(e) + inf ip*(a) = 0. 

eev v aev 1 - 

We conclude this section extending these concept to integral functionals. Note 
that if ip : M. m x Jl — >■ R satisfies assumptions (1.2), then for all u G L P (Q) the map 
u> \— ?■ ip(u(u>),uj) is measurable and 



$(«):= / ip(u(co), co) dfi < oo. 
Jn 

We can now connect integral functionals and conjugate functions with the following 
theorem. 
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Theorem 2.4 (Rockafellar [25]). Let ip : IR m x Q — >■ R satisfy assumptions (1.2), 

$*(z) = / Lp*(z(u),u)d[i. 



n 



3 Proof of Theorem 1.1 

Proof. The proof is divided into three steps. 

Step I. The auxiliary problem. Let £ G R m be fixed, we look for functions 
v G L p (tt), z G L p \Q) such that 

vev; ot , zel/jfl) with z ed^ + v). (3.14) 

By Lemma 2.6-(iii) and (iv), (3.14) is equivalent to 

veV? ot , zeiC(n) with t + v ed<p*{z). (3.15) 

Define the functionals 

F:L P {Q)^R, F(u):= [ </>(£ + u(u),u)dii, (3.16) 

Jq 



G : V (Q.) -* R, G(y) := / (/>*{y(u),u) - £ • y(w) d/x. (3.17) 

Jn 

We claim that (3.14) (or (3.15)) is equivalent to the minimization problem: 

v solves: min F(u) (3.18a) 



«6V£ ot 



and 



z solves: min G(y). (3.18b) 

Let us show that (3.14) implies (3.18a). For / G L\Q), let / := / - E(f) be the 
fluctuation of /, so that 

E(f) = 0, f = f + E(f). 
By definition of subdifferential, z G <90(£ + v) if and only if 

z-(\-^ + v)) + ^ + v) <0(A) VAel™ 
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Then, for all u G Lp Qt (Q) such that E{u) = £, we have, /i-a.e. in Q, 

z-(Z + u-(£ + v)) + <j>(Z + v)< 0(£ + u). 
By Weyl's Lemma 2.4, it holds 



z(uj) • u{uj) dpi — — / z(ui) • v(u>) d/j,, 
n Jn 

and we conclude 

0(£ + v(u>), u) dji < I <j>(£ + u(oS),u) d/j, 



for all u G Lp 0t (f2) such that E{u) = 0, i.e., for all u G Vp 0t . The proof that (3.15) 
implies (3.18b) follows in the same way. 

Let us prove the converse implication. By definition of F (3.16) and definition 
of conjugate function (2.12), for all y G L p (fl) it holds 

F*(y)= sup {(y,u)-F(u)} 
ueL p (n) 



sup < / y ■ u — 0(£ + u) d\x 

ueL p (n) {Jn 



sup < y ■ (v - £) - 0(w) d/x 

sup <^ (y, u) - / 0(u) dju 
«eL p (n) I Jn 



Th.2A 



'(n) (. Jn ) Jn 



4>*{y) -y-£dfi = G(y). 

n 



Let now £ G M m and let (v,z) be a solution of (3.18). Since by Weyl's theorem 
(Vpot) -1 = Lg 0l (Q), by Theorem 2.3 we obtain 



/ T 1 /? 2 ^ 

/ 0(C + t ') = mm e(u) = ~~ mm F*{y) 

Jn ' «e V Dot «ez/,(m 



-G(*) = U-z-<P{z)= f{^ + v)-z-<t>*{z). 
Jn Jn 

(3.19) 



Since by Lemma 2.6-(ii) 

0(£ + f) > (£ + v) • -2 — <f>*(z) for /i-a.e. cu G Jl, 
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by (3.19) we conclude that 

0(£ + v) + <f>*(z) — (£ + v) • z for /i-a.e. u G f2, 

and therefore, by (2.6)-(v), that (v,z) solves (3.14). 

We turn now to the existence of a solution to the minimization problems (3.18a)- 
(3.18b). We note that 

• by Lemma 2.3 the subspaces L P ol (Q) and Vp 0t are closed with respect to the 
weak convergence in L P (Q). 

• 0,0* are convex functions, with growth (1.2c) and (2.13) respectively, which 
implies that the functionals F and G defined in (3.16) and (3.17) are coercive 
and lower-semicontinuous with respect to the weak convergence in L P (Q). 

Then, any minimizing sequence for F (G) is weakly compact in L P (Q) (resp. 
LP (Q)) and, owing to lower-semicontinuity, any weak limit point is a solution 
of (3.18a) (resp. (3.18b)). 

Remark 3. If and 0* are strictly convex, then the solution of (3.18) is unique. 
Under the hypothesis of Theorem 1.1 problem (3.18) might have more than one 
solution. 

Step II. The homogenized functional 0o- Let us define 0o : M m — > K 



O (O:= inf / 0(£ + u{u),u)dn. (3.20) 

«ev£ ot J n 

We claim that 

for all £ G M m , every solution (v, z) of problem (3.14) satisfies 

E{z) e dMO- (3 ' 21) 

The proof is a direct translation to functions in L P (Q) of the proof of Theorem 3.2 

in [34]. 

Proof of (3.21). Recall that E(f) denotes the expected value f n f d\i of a random 
function /. Define ip : R m -> R 

MC)--=^yj*(z(u),co)d^ : zel£ A (n),E(z)={\. (3.22) 
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It is not obvious, a priori, whether ipo = 0q or not. It will become evident by the 
end of the proof, but at the moment we have to treat this two functions as separate 
objects. 

Owing to Fenchel's inequality (Lemma 2.6-(ii)), VA G W n , Vw G Vp 0t , V s G 



X ■ E(s) = (X + u(uj)) ■ s(oo) dfi < / 0(A + u(u),u) + 4>*(s(ui),uj) dfj,. 
Jn Jn 

Passing to the infimum over u G Vp 0t and recalling definitions (3.20) and (3.22) we 
get 

A-£(s)<0o(A)+Vo(£(s)), (3.23) 

and passing to the supremum on A G W 71 we read 

<i>*(E(s)) = sup {A • E(s) - O (A)} < ME(s)). (3.24) 

AeM m 

Let now v,z be the solutions of the auxiliary problem 

"£V; t) ?6Li(fi) with ^G90(e + t;), (3.25) 

for C, G M m . By Fenchel's equality (Lemma 2.6-(v)) 

£-E(z)= I \i + v(u)) ■ z(u) dfi = [ (f>(t + v(ou),ou) + (f)*(z(ou),uj)dii, 
Jn Jn 

and therefore, by definitions (3.20) and (3.22), we get 

£-E(z)>MS)+ME(z)). (3.26) 

Combining (3.23) and (3.26) we obtain that for all z which solve (3.25) 

ME(z))=<l>* (E(z)). (3.27) 

Combining (3.24), (3.26), and (3.27), we obtain that for all f G R m , for all (v,z) 
satisfying (3.25), 

^■E(z)=MO + ^o(E(z)), 

and therefore, by Fenchel's equality, that 

E(z) G dMO- 
This concludes the proof of (3.21) and Step II. □ 

Additionally, we note that since for all A G lR m there exist z G L^ ol (Q), u G 
Lp 0t (fi) such that 

E(z) = A, z G d(f)(u), 

equality (3.27) is valid for all A G M. m , and therefore ipo = 4>q. 
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Step III. Convergence. We recall that an operator A : D(A) C X — ¥ V(X), 
defined on a Hilbert space X, is said to be monotone if 

(y2-yi,x 2 -x 1 )>0 VxiED(A), V yi eA(xi). 

A monotone operator is said to be maximal monotone if its graph is maximal 
among all monotone sets, in the sense of set inclusions. We are going to use the 
following two fundamental properties. 

• The subdifferential of a convex function : IR m — > R is a maximal monotone 
operator. 

• Let x, y G M m , if d<j) is maximal monotone, then 

JVxGM™, 3ye d(f)(x), (y-y,x-x)>o\ =4> y G d(j>{x). (3.28) 

Let now £,v,z be as in (3.25), we define the oscillating test functions 

v v (x) := x H- v (T x/v u) , z v (x) := x H- z (T x/v uj) , 

<j> v (;x) = <f)(-,T x/v uj), 
which satisfy 

v v eL p pot (R m ;R m ), (3.29a) 

^GL^(M m ;M m ), (3.29b) 

Zr,ed^ + v v ). (3.29c) 

We can find a set N C Q, with //(iV) = 0, such that conditions (1.2b) and (1.2c) 
are satisfied for all u G Q\N. Fix u G Q\N and consider now the problem 

-diva" = P 

o- r > G <90„(W) 
m" = 

As in the proof of Step II above, (a^^u 71 ) is a solution of (3.30)-(3.31) if and only 
if u v realizes 



weakly in Q, 


(3.30) 


a.e. in Q, 


(3.31) 


a.e. on dQ. 


(3.32) 



mm 



(f>r,(Vu(x),x) - P{x)u{x) dx, u G Wq' p (Q) 1 (3.33) 

> J 
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and a v realizes 



mm 



| f (f>* v (s(x),x)dx, seL p '(Q;R m ), -div s = f 1. (3.34) 



By direct methods of calculus of variations, as in the proof of Step II, for all 77 > 
there exists a (not necessarily unique) solution {vP^a 11 ). Denoting by p (, ) p > the 
duality between W ' P (Q) and W~ ,p (Q), by Lemma 2.6 -(v), (1.2c), and Lemma 
2.5 

minKcJ (||V«1*, W) + lkl£, (giRm) ) " CblQl 

< f 0„(Vt/ ; ) + 0;(a") (fa = /" Vu" >a«dx = p (u , /"),,. 



Since (f v ) is a converging sequence, we can assume that ll/^llw-i.pVn) * s bounded by 
a constant C/, and by Poincare's lemma we can find a constant C > 0, depending 
on Q, C/, Co, ci, ci, but independent of 77, such that 

By compactness, for all u G f2\iV we can find a subsequence (77^) and a limit point 
(u*(w),<T*(a;)) = (u*,a*) G Wq' p (Q) x L p ' (Q;R m ) such that 

u"(-,w) -^ u* weakly in W^ P (Q) (3.35) 

and 

^(-,00)^(7* weakly in L P '{Q). (3.36) 

Testing equation (3.30) with g G VF ' P (Q) we have 

/ <r«.Vgdx= P (g,n p >, 

JQ 

and passing to the limit as 77 — > we find 

/ a* -Vgdx= p {g, f) p/ , 
JQ 



that is, 



It remains to prove that 



-div<r* = /. 

a* G <90 o (Vw*) 
19 



By monotonicity of d<f) v and equations (3.31) and (3.29c), \/i] > 0, V£ G M. m , for 
a.e. x G Q, and /i-a.e. cu G fi it holds 

(a" - z v ) ■ (Vu" - (£ + v„)) > 0. (3.37) 

By Birkhoff 's ergodicity Theorem 2.2 and the definition of v v and ^ 

^ -^ E(z) weakly in L P '(Q) and t,+v v -±{, weakly in L P (Q). (3.38) 

By (3.30) and (3.29b) 

{divo"''} , {divz T; } are compact sets in W~ ,p (Q), 

and by (3.29a) curlf^ = 0. Therefore the hypothesis of the div-curl Lemma 2.2 are 
satisfied, and by convergences (3.35), (3.36), and (3.38) we can pass to the limit 
as i] — > 0, in (3.37), obtaining 

(a* - E(z)) ■ (Vm* - > for a.e. xeQ. 



m 



Owing to Step II, we know that E(z) G <9</>o(£), and by the arbitrariety of £ G 
and the maximality property (3.28) we obtain that a* G (90o(Vw*) for a.e. x G 
This concludes the proof of Theorem 1.1. □ 



References 

[1] Y. Abddaimi, G. Michaille, and C. Licht. Stochastic homogenization for an 
integral functional of a quasiconvex function with linear growth. Asymptot. 
Anal, 15(2):183-202, 1997. 

[2] I. Babuska. Homogenization and its application. Mathematical and computa- 
tional problems. In Numerical solution of partial differential equations, III 
(Proc. Third Sympos. (SYNSPADE), Univ. Maryland, College Park, Md., 
1975), pages 89-116. Academic Press, New York, 1976. 

[3] A. Bensoussan, J.-L. Lions, and G. Papanicolaou. Asymptotic analysis for 
periodic structures, volume 5 of Studies in Mathematics and its Applications. 
North- Holland Publishing Co., Amsterdam, 1978. 

[4] A. Bourgeat, A. Mikelic, and S. Wright. Stochastic two-scale convergence in 
the mean and applications. J. Reine Angew. Math., 456:19-51, 1994. 

20 



[5] L. A. Caffarelli and P. E. Souganidis. Rates of convergence for the homog- 
enization of fully nonlinear uniformly elliptic pde in random media. Invent. 
Math., 180(2):301-360, 2010. 

[6] L. A. Caffarelli, P. E. Souganidis, and L. Wang. Homogenization of fully 
nonlinear, uniformly elliptic and parabolic partial differential equations in 
stationary ergodic media. Comm. Pure Appl. Math., 58(3):319-361, 2005. 

[7] G. Dal Maso and L. Modica. Nonlinear stochastic homogenization. Ann. Mat. 
Pura Appl. (4), 144:347-389, 1986. 

[8] G. Dal Maso and L. Modica. Nonlinear stochastic homogenization and ergodic 
theory. J. Reme Angew. Math., 368:28-42, 1986. 

[9] A. Damlamian, N. Meunier, and J. Van Schaftingen. Periodic homogenization 
for convex functionals using Mosco convergence. Ric. Mat., 57(2):209-249, 
2008. 

[10] E. De Giorgi and S. Spagnolo. Sulla convergenza degli integrali dell'energia 
per operatori ellittici del secondo ordine. Boll. Un. Mat. Hal. (4), 8:391-411, 
1973. 

[11] J. L. Doob. Stochastic processes. Wiley Classics Library. John Wiley & Sons 
Inc., New York, 1990. Reprint of the 1953 original, A Wiley-Interscience 
Publication. 

[12] I. Ekeland and R. Temam. Analyse convexe et problemes variationnels. Dunod, 
1974. Collection Etudes Mathematiques. 

[13] A. Gloria and F. Otto. An optimal variance estimate in stochastic homoge- 
nization of discrete elliptic equations. To appear on Annals of Probability. 

[14] V. V. Jikov, S. M. Kozlov, and O. A. Olemik. Homogenization of differential 
operators and integral functionals. Springer- Verlag, Berlin, 1994. Translated 
from the Russian by G. A. Yosifian [G. A. Iosifyan]. 

[15] S. M. Kozlov. The averaging of random operators. Mat. Sb. (N.S.), 
109(151)(2):188-202, 327, 1979. 

[16] P.-L. Lions and P. Souganidis. Homogenization for stochastic viscous 
Hamilton- Jacobi equations. To appear on Comm. Partial Differential Equa- 
tions. 



21 



[17] K. Messaoudi and G. Michaille. Stochastic homogenization of nonconvex in- 
tegral functional. Duality in the convex case. Sem. Anal. Convexe, 21:Exp. 
No. 14, 32, 1991. 

[18] K. Messaoudi and G. Michaille. Stochastic homogenization of nonconvex in- 
tegral functionals. RAIRO Model. Math. Anal. Numer., 28(3):329-356, 1994. 

[19] F. Murat. Compacite par compensation. Ann. Scuola Norm. Sup. Pisa CI. 
Sci. (4), 5(3):489-507, 1978. 

[20] F. Murat. Compacite par compensation: condition necessaire et suffisante de 
continuity faible sous une hypothese de rang constant. Ann. Scuola Norm. 
Sup. Pisa CI. Sci. (4), 8(1):69-102, 1981. 

[21] G. Nguetseng. A general convergence result for a functional related to the 
theory of homogenization. SI AM J. Math. Anal, 20(3):608-623, 1989. 

[22] G. C. Papanicolaou and S. R. S. Varadhan. Boundary value problems with 
rapidly oscillating random coefficients. In Random fields, Vol. I, II (Eszter- 
gom, 1979), volume 27 of Colloq. Math. Soc. Jdnos Bolyai, pages 835-873. 
North-Holland, Amsterdam, 1981. 

[23] G. C. Papanicolaou and S. R. S. Varadhan. Diffusions with random coef- 
ficients. In Statistics and probability: essays in honor of C R. Rao, pages 
547-552. North-Holland, Amsterdam, 1982. 

[24] F. Peter and H. Weyl. Die Vollstandigkeit der primitiven Darstellungen einer 
geschlossenen kontinuierlichen Gruppe. Math. Ann., 97(l):737-755, 1927. 

[25] R. T. Rockafellar. Integrals which are convex functionals. Pacific J. Math., 
24:525-539, 1968. 

[26] E. Sanchez-Palencia. Nonhomogeneous media and vibration theory, volume 
127 of Lecture Notes in Physics. Springer- Verlag, Berlin, 1980. 

[27] B. Schweizer. Homogenization of the Prager model in one-dimensional plas- 
ticity. Contin. Mech. Thermodyn., 20(8):459-477, 2009. 

[28] B. Schweizer and M. Veneroni. Periodic homogenization of the Prandtl-Reuss 
model with hardening. J. Multiscale Modeling, 2(l-2):69-106, 2010. 

[29] N. Svanstedt. Multiscale stochastic homogenization of monotone operators. 
Netw. Heterog. Media, 2(1):181-192 (electronic), 2007. 

22 



[30] N. Svanstedt. Multiscale stochastic homogenization of convection-diffusion 
equations. Appl. Math., 53(2):143-155, 2008. 

[31] L. Tartar. Problemes de controle des coefficients dans des equations aux 
derivees partielles. In Control theory, numerical methods and computer sys- 
tems modelling (Internat. Sympos., IRIA LABORIA, Rocquencourt, 1974), 
pages 420-426. Lecture Notes in Econom. and Math. Systems, Vol. 107. 
Springer, Berlin, 1975. 

[32] L. Tartar. Cours Peccot au College de France. Partially written by F. Murat 
in Seminaire d'Analyse Fonctionelle et Numerique de l'Universite d'Alger, 
unpublished, 1977. 

[33] L. Tartar. Compensated compactness and applications to partial differential 
equations. In Nonlinear analysis and mechanics: Heriot-Watt Symposium, 
Vol. IV, volume 39 of Res. Notes in Math., pages 136-212. Pitman, Boston, 
Mass., 1979. 

[34] A. Visintin. Scale-integration and scale-disintegration in nonlinear homoge- 
nization. Cole. Var. Partial Differential Equations, 36(4):565-590, 2009. 



23 



